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THE REMARKABLE 7-ray ACTIVITY IN THE GRAVITATIONALLY LENSED BLAZAR PKS 1830-211 
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ABSTRACT 

We report the extraordinary 7-ray activity (E > 100 MeV) of the gravitationally lensed blazar PKS 
1830-211 (z = 2.507) detected by AGILE between October and November 2010. The source experi- 
enced on October 14 a flux increase of a factor of ~ 12 with respect to its average value and kept 
brightest at this flux level (~ 500 x 10~ 8 photons cm~ 2 s~ 4 ) for about 4 days. The 1-month 7-ray 
light curve across the flare showed a mean flux ¥(E > 100 MeV)= 200 x 10~ 8 photons cm -2 s _1 , 
which resulted in an enhancement by a factor of 4 with respect to the average value. Following the 
7-ray flare, the source was observed in NIR-Optical energy bands at the Cerro Tololo Inter- American 
Observatory and in X-rays by Swift/XRT and INTEGRAL/IBIS. The main result of these multifre- 
quency observations is that the large variability observed in 7-rays has not a significant counterpart 
at lower frequencies: no variation greater than a factor of ~ 1.5 resulted in NIR and X-ray energy 
bands. PKS 1830-211 is then a good "7-ray only flaring" blazar showing substantial variability only 
above 10-100 MeV. We discuss the theoretical implications of our findings. 

Subject headings: galaxies: active; galaxies: jets; quasars: general; quasars individual (PKS 1830-211); 
radiation mechanisms: non-thermal 



1. INTRODUCTION 

PKS 1830-211 is a high redshift blazar (z = 2.507, Lid- 
man et al. 1999) gravitationally lensed by a spiral galaxy 
at z = 0.886 (Wiklind & Combes 1996), as showed by the 
two radio lobes located 1" apart from each other (Lovell 
et al. 1998). The lensed counterparts were observed also 
in near infrared (NIR) and optical energy bands by Hub- 
ble Space Telescope and Gemini observatories (Courbin 
et al. 2002). In X-rays the source was observed by 
both XMM-Newton and Chandra allowing to study in 
details the complicated soft X-ray behaviour, probably 
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due to absorption (Dai et al. 2008). Moreover, the source 
presents a soft 7-ray energy spectrum (L = 2.56) as de- 
tected by EGRET (Hartman et al. 1999) and recently 
confirmed by Fermi LAT (Abdo et al. 2010). 

The modelling of the SED of this blazar (De Rosa et al. 
2005) led to the tentative classification of PKS 1830-211 
as a MeV blazar (the inverse Compton peak lies below 
100 MeV, Sikora et al. 2002), with the caveat that the 
whole dataset was not simultaneous. In particular, the 
SED was interpreted in the standard framework of a one- 
zone leptonic model where the seed photons responsible 
for inverse Compton peak should be originated within 
the torus surrounding the central engine. 

The source has been extensively monitored in hard X- 
rays by INTEGRAL/IBIS on 7-year long period, dur- 
ing which no significant variability has appeared on both 
short (daily) and longer (monthly) timescales (Zhang et 
al. 2008). 

Since its launch, AGILE (Tavani et al. 2009) detected 
a 7-ray flare (F(E > 100 MeV)= (160 ± 50) x 10" 8 ph 
cm" 2 s- 1 ) of PKS 1830-211 in October 2009 (Striani et 
al. 2009) which resulted in an enhancement by a factor of 
3 with respect to the 1-week average flux detected before 
the flare (F(E > 100 MeV)= 60 x 10~ 8 photons cm" 2 

One year later, Fermi LAT detected a significant en- 
hancement of the 7-ray emission F(E > 100 MeV)= 
(520 ± 110) x 10~ 8 photons cm" 2 s" 1 and F(E > 100 
MeV)= (1400 ±500) x 10~ 8 photons cm~ 2 s" 1 on 1-day 
and 6-hr timescales, respectively (Ciprini et al. 2010), 
resulting in an increase of a factor of about 12 and 35 
with respect to the average value reported by Abdo et 
al. 2010. Moreover, AGILE detected a prolonged 7-ray 
activity two days after the flare reported by Fermi (Don- 
narumma et al. 2010). 
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In this letter, we present the results obtained by means 
of a multifrcquency monitoring of PKS 1830-211 from 
NIR to 7-ray energy bands. In detail, we describe the 
multi-A data analysis and a possible interpretation of 
the SED of this object to reconcile the average 7-ray 
activity with the 7-ray enhancement. We adopted a A- 
CDM cosmology with the following parameters h = 0.71, 
n m = 0.27, fl A = 0.73. 

2. THE MULTIWAVELENGTH CAMPAIGN 
2.1. AGILE Observations 

AGILE/GRID data were analyzed using the Build-20 
software. Well-reconstructed 7-ray events were selected 
using the FM3.119 hltcr. All the events collected during 
the passage in the South- Atlantic Anomaly were rejected. 
We filtered out the Earth-albedo, rejecting photons com- 
ing from a circular region of radius 80 degrees and cen- 
tered on the Earth. A 1-month long 7-ray light curve (2- 
day time-bin, see Fig. 1) across the 7-ray flare reported 
by Fermi/LAT has been produced between 2010-10-08 
00:00UT and 2010-11-08 00:00UT (MJD 55477-55509) 
using the standard AGILE Maximum-Likelihood proce- 
dure. We obtain a mean flux on the whole integration 
time of F(E > 100 MeV)=(200 ± 29) x 10~ 8 photons 
cm~ 2 s _1 which is a factor of ~ 4 higher than the flux re- 
ported in the first Fermi catalog. We note that this 7-ray 
enhancement represents a 5-a excess with respect to its 
value in steady state as detected by AGILE (F(E > 100 
MeV)= (32 ±8) x 10~ 8 photons cm" 2 s"\ Verrecchia et 
al. in preparation). 

We extracted the 7-ray spectrum between lOOMeV — 
lGeV, discarding the energy channels above 1 GeV due 
to the poor statistics. The AGILE spectrum is well fitted 
by a power-law with a photon index of 2.05±0.17 (all the 
errors reported in this paper are at la, unless otherwise 
stated). 

As shown in Fig.l, the source remained at a flux level 
~ 200 x 10 -8 photons cm~ 2 s _1 on a 1-month time 
scale. A clear enhancement emerged on October 14 
(MJD 55483) that lasted about 4 days. The 2-day binned 
AGILE data showed that the source reached its maxi- 
mum of F(E > 100 MeV)~ (500 ± 130) x 10~ 8 photons 
cm" 2 s- 1 on October 15 (MJD 55484.75 ± 0.7); the esti- 
mate of this epoch accounts for the bias introduced by the 
choice of the starting time (T s t or t) and binning adopted 
to obtain the light curve. This was done by operating 
4 translations of T start (by 0-, 0.5-, 1-, 1.5-day) in the 
light curve. 

We extracted the 7-ray spectrum also during the 4-day 
flare, obtaining a power-law photon index V = 2.40±0.28. 
Because of the poor statistics on this short duration, 
it was not possible to infer any hardening/softening oc- 
curred during this flare, being the photon indices com- 
patible within 1-er. In section 3, we will focus on the 
1-month 7-ray enhancement and on the 4-day flare de- 
tected between October 14 and 18 (horizontal lines in 
Fig. 1). 

2.2. X-ray observations 

2.2.1. The INTEGRAL observations 

During the 7-ray flare of PKS 1830-211 in October 
2010, INTEGRAL' (Ubertini et al. 2003) was monitor- 
ing the Galactic bulge region (since INTEGRAL AO-3, 



Kuulkers et al. (2007) have been observing this region 
regularly during all the visibility periods) . The monitor- 
ing across the 7-ray flare (MJD 55479 - 55491, 2010 Oc- 
tober 10-22) did not allow a detection of PKS 1830-211, 
then an upper limit of^3xl0~ n erg cm -2 s" 1 (it ac- 
counts for the 15 degrees off-axis position of the source) 
was obtained on a net exposure of 200 ksec (October 
14-18). Since the hard X-ray average flux (20-40 keV) 
is 2.6 mCrab (1.95 x lO" 11 erg cur 2 s"\ Bird et al. 
2010), we may exclude variations greater than a factor 
of 1.5 in this energy band. Therefore, in this paper we 
adopted the INTEGRAL/IBIS data collected in the 4th 
IBIS survey (Bird et al. 2010) consisting of 7355 Science 
Windows (lasting about 2000 s each), performed from 
the beginning of the Mission in November 2002 up to 
April 2008 and including all available public and Core 
Programme data. The total on-source exposure is 4258 
ks. IBIS/ISGRI images for each available pointing were 
generated in various energy bands using the ISDC off- 
line scientific analysis software OSA (Goldwurm et al. 

2003) version 7.0. Then all images have been mosaicked 
to create significance map at revolution level (each rev- 
olution lasting about 3 days) and for the all available 
pointings. Count rates at the position of the source were 
extracted from individual images in order to provide light 
curves in various energy bands; from these light curves, 
we extracted and combined the average fluxes in order 
to produce an average source spectrum (see Bird et al. 
2010 for details). The derived hard X-ray spectrum is 
well fitted by a power-law with T = 1.56 ± 0.14, provid- 
ing a flux of 1.4 x 10~ n erg cm~ 2 s _1 and 2.7 x 10~ n 
erg cm~ 2 s _1 in 20-40 keV and 40-100 keV, respectively. 
It is worth noting that this spectrum represents the best 
coverage of PKS 1830-211 in the hard X-rays. 

2.2.2. The Swift observations 

Following the 7-ray flare, Swift/XRT (Gehrels et al. 

2004) pointed at PKS 1830-211 twelve times, between 
2010 October 15-27 (MJD 55484-55496) with net ex- 
posures ranging between 1000 and 4000 s. The XRT 
level 1 data were processed with the standard procedures 
(XRTPIPELINE v. . 12 . 2). We selected photon counting 
(PC) mode data with the standard 0-12 grade selection, 
and source events were extracted in a circular region of 
aperture ~ 40", the background was estimated in differ- 
ent same sized circular regions far from the source. We 
created response matrices through the xrtmkarf task. 
As for spectra with very few counts, we applied the Cash 
statistics and we rebinned the remaining spectra in or- 
der to have at least 30 counts per bin, excluding from 
the fit the channels with energies below 0.3 keV and 
above 10 keV. We fitted each spectrum, with a contin- 
uum power-law absorbed both with a Galactic column 
density (A^^f al =2.6 x 10 21 cur 2 Stark et al. 1992) and an 
additional absorber located at a redshift z=0.886. Due 
to the short exposure for each observation, the fitting 
procedure resulted in a degeneracy between and T. 
The values of we found combining all the spectra 
together are in full agreement with the previous mea- 
surement obtained with Chandra observations, therefore 
we decided to fix 7V| to this value, i.e. 1.94xl0 22 cm -2 . 
All data sets are well reproduced by this model, with 
X 2 /dof ranging between 0.9 and 1.1. The best fit values 
of the photon index range between r=1.0±0.3-1.4±0.2 
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with absorbed fluxes in 0.3-10 keV between 0.91q'j and 



1.2lSixlO- u 



erg cm 2 s , all consistent within 1-a 



given the large uncertainties. On 2010 October 15, the 
day overlapping the gamma-ray flare, the absorbed flux 
in '0.3-10 keV is l.l±g:|xl0 _ " erg cm" 2 s" 1 at 90% 
confidence level. We note that these fluxes are well con- 
sistent with the ones reported by De Rosa et al. (2005), 
ruling out any significant variation occurred during the 
high 7-ray activity (Fig. 1). On the time interval of 
the gamma-ray enhancement (15-18 oct. 2010), the full 
band detected count rates are between 0.095j^oo6 an( ^ 



0.148_ 009 counts/s, showing a variation amplitude that 
is less than a factor of 1.6. 

2.3. Optical monitoring: SMARTs Consortium 

Optical and infrared data on PKS 1830-211 were ob- 
tained using the 1.3m telescope at the Cerro Tololo Inter- 
American Observatory under the Small and Moderate 
Aperture Research Telescope System (SMARTS) pro- 
gram. We obtained simultaneous data in the optical 
R and infrared J bands using the ANDICAM instru- 
ment, allowing the acquisition of data from 0.4 to 2.2 
microns. Our dataset consists of 8 i?-band and 17 J- 
band images, taken between 2010 October 16 and 27 
(MJD 55485-55496). Images were taken approximately 
every three nights, with exposure times of 400 to 600 
seconds in R, and 90 to 600 seconds in J. Both the op- 
tical and IR data were flat-fielded, bias-subtracted, and 
over-scan corrected using standard IRAF packages. Sin- 
gle exposures were not sufficient to detect the object in 
either waveband, nor were summations of the images us- 
ing standard IRAF packages. We estimate a limiting 
magnitude of 21.5 in R and 17.6 in J, based on the com- 
bined images and the magnitudes of field stars acquired 
from the 2MASS. The observations performed during the 
4-day gamma-ray flare resulted in a limiting magnitude 
17.6 and 19 in J and R, respectively. The J magni- 
tude limit during the flare period remained unchanged 
from the longer summation due the sky limit already 
reached by the shorter exposure time. We note that our 
limit magnitude could be contaminated by the star SI 
at 1" from PKS 1830-211 (Meylan et al. 2005) unre- 
solved by SMARTs. Neverthless, no enhancement at or 
above two standard deviations of the background (sky) 
variance was detected within 2" of the position of PKS 
1830-211. Therefore, by comparing our limiting magni- 
tude in J with the photometric data of PKS 1830-211 
(J = 18.04), we exclude variation > 1.5 in J-band. 

3. DISCUSSION 

Here we discuss the SED modelling of the blazar PKS 
1830-211 during the high 7-ray state observed between 
October-November 2010. Multifrequency observations 
were activated to follow the 7-ray flare covering NIR- 
Optical (SMARTs), soft and hard X-ray energy bands 
(Swift/XKT and INTEGRAL/ISGRI, respectively). The 
main results obtained during this multifrequency cam- 
paign are that the high and unusual activity recorded in 
7-rays seems to have no significant counterpart at lower 
frequencies. In detail, the simultaneous NIR-optical, soft 
and hard X-ray emissions of this source did not follow the 
significant changes observed in 7-rays. As for the NIR 
and optical emissions, we derived only upper limits for 



this source that is indeed vary faint especially in opti- 
cal energy bands (Courbin et al. 2002). Neverthless, 
by comparing the derived upper limits in J-band (red 
downward arrow in Fig. 2) with the NIR data reported 
by Courbin et al. 2002 and Meylan et al. 2005, we can 
exclude variation in the thermal + non-thermal compo- 
nents greater than a factor of ~ 1.5. The soft and hard 
X-ray emissions also exclude variation greater than a fac- 
tor of 1.6 and 1.5, respectively. It is worth noting that the 
observed variations of the SED rule out the hypothesis 
that the 7-ray emission was connected to macrolensing, 
since its effects would be energy- independent. On the 
other hand, the chromaticism of the SED variability may 
suggest that microlensing from stars in the lensing galaxy 
may cause the observed gamma-ray variability (see e.g. 
Torres et al. 2003). In order to investigate this possibil- 
ity, we simulated the transition of a source with different 
dimensions (Jovanovic et al. 2008) across the microlens- 
ing magnification pattern for the lens system (taking con- 
vergence k = 0.158 and shear 7 = 0.096 as given in Winn 
et al. 2002). By assuming a gamma-ray emitting region 
of 10 15 cm (observer frame), we found an amplification by 
a factor of tens on time scale of years. Even worst is the 
case of the gamma-ray flare observed by Fermi on 6-hr 
time scale, which is a factor of 14 higher than the steady 
state, requiring high amplification in very short time. 
Shorter time-scales of microlensing can be obtained with 
very compact source (order of 10 13 cm), but we cannot 
further reduce the gamma-ray emitting region otherwise 
it would result in an increase in gamma-ray opacity due 
to pair production. Therefore, we found that the ob- 
served variability is more likely intrinsic to this blazar. 

We complemented the low frequency data with the 
non-simultaneous Planck observations (30, 100, 217 
GHz) taken from the Early Release Compact Source Cat- 
alogue (ERCSC) extracted from the first all-sky survey of 
Planck 17 . In Fig. 2 we report all the observed multifre- 
quency data and the SED modelling taking into account 
the acromatic magnification due to the gravitational lens- 
ing, assumed to be one order of magnitude (De Rosa et 
al. 2005). 

We model three 7-ray states of the source. The first 
one is the average state reported by De Rosa et al. 2005 
and consistent with the average 7-ray state reported in 
Abdo et al. 2010 (~ 40 x 10~ 8 ph cm" 2 s" 1 ) (see black 
points in Fig. 2); the second is the 1-month enhance- 
ment showing a mean flux of ~ 200 x 10~ 8 ph cm~ 2 
s _1 (blue points in Fig. 2); the third is the peak of the 
7-ray emission (October 14-18; hereafter, "flare") which 
is a factor of ~ 2 — 3 higher than the 1-month average 
flux (red points in Fig. 2). We note that the multifre- 
quency behavior of the source, i.e. the lack of correlated 
variability between the low (NIR-optical bands, X-rays) 
and high energy (7-rays) portions of the SEDs, disfavors 
the one-zone leptonic model for this event. We account 
for this evidence, assuming two relativistic electron pop- 
ulations with the energy density described by a broken 
power law in the Lorentz factor domain (7) with break 
equal to jb- 



3(7) 



(1) 



(l/lb) ai + {l/lb) ah 

1 http : //www. sciops . esa. int/ index. php?project=planck&page=Planck_Lej 
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where a; and ah are the pre- and post-break electron 
distribution spectral indices, respectively. We assumed 
that the two electron populations contained a random 
magnetic field B. In Table 1 we summarize the values of 
the model parameters for the two components. 

We estimated the contributions of both electron pop- 
ulations to sychrotron emission and to synchrotron self 
Compton (SSC) plus inverse Compton from external seed 
photons (EC) originating in a dusty torus, in the accre- 
tion disk and in the Broad Line Region (BLR) (see Sikora 
et al. 2002, Vittorini ct al. 2009). In detail, we assume: 
1) the emission from a torus with temperature T = 10 
K distant 1.26 x 10 19 cm from the central engine 2); the 
direct emission from the disk with T = 1.5 x 10 4 K and 
luminosity Ldisk — 6 x 10 45 erg s _1 distant 5 x 10 16 cm 
from the 2nd population; 3) the emission from the BLR 
distant Rblr = 5 x 10 17 cm from the central engine and 
reprocessing 10% of the disk luminosity. 

In this scenario, the average 7-ray state is produced 
through EC on dusty torus, disk and BLR seed photons 
mainly by the first population, resulting in a jet power 
Pj et = 10 46 erg s^ 1 (black line in Fig. 2). The 1-month 
7-ray enhancement, a factor of 4 higher than the aver- 
age flux, has required the additional contribution of the 
second electron population, which is characterized by a 
smaller size, higher electrons density and higher 7^ with 
respect to the 1st one (Table 1). On the basis of these as- 
sumptions, the 7-ray enhancement (blue points in Fig. 2) 
is then produced by EC of the 2nd population on the 
photon density field while moving inside the BLR with 
T = 15, transporting a jet power Pj et = 3 x 10 46 erg s _1 . 
This modelling is in agreement with the source fading 
at its average flux level when the blob is moving outside 
the BLR (At ~ Rblr{1 + z)/c8 > 30 days). Follow- 
ing this description, the synchrotron and IC emissions 
are dominated by the 1st population, thus accounting 
for the moderate variation inferred by the NIR-optical 
and X-ray observations. On the other hand the "flare" , 
which is a factor of 12 greater than the average flux and 
a factor of ~ 2 — 3 greater than the 1-month activity, re- 
quires a "local" enhancement of the photon density field 
by a factor of 3 (red line in Fig. 2), likely due to a blob- 
cloud interaction (Araudo et al. 2010). In this case we 
derived a jet power of Pj l ^ re = 6 x 10 46 erg s" 1 . The 
lack of correlated variabilities in optical, X- and 7-rays 
during the flare prevented us to explain this event as as- 
sociated with changes of the electron population in terms 
of injection/acceleration mechanisms. 
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We underline that the 7-ray behavior of PKS 1830-211 
recorded by AGILE and Fermi/LAT between October- 
November 2010 is rare as probed by the lack of similar 
variability since 2007. Similarly, no hard X-ray variabil- 
ity has been observed during the 7-year monitoring by 
INTEGRAL/IBIS. In conclusion, we can attribute the 
properties of PKS 1830-211 to the ones of "7-ray only 
flaring" blazars showing substantial variations especially 
in the 7-ray energy range above 10-100 McV. These evi- 
dences strongly support our interpretation of a "steady" 
electron population, filling the jet below the BLR, that 
is responsible for the average 7-ray emission recorded 
by both EGRET and Fermi/LAT, more than ten years 
apart. 

It is worth noting that the inspection of the 7-ray light 
curve on a longer period across the October-November 
enhancement, did not show any evidence of the time de- 
lay between the emissions of the two lensed images A and 
B as measured in the radio maps (26^4 days, Lovell et al. 
1998). We can assume the lack of the delay between A 
and B only if the flux ratio of the 2 components is <~ 1 . If 
it was lower than 1 , we are prevented to draw any conclu- 
sion because the emission would be occurred below the 
AGILE sensitivity level. We also note that microlensing 
could play a role in explaining the possible lack of echo 
in the gamma-ray light curve. In fact, in some lensed 
quasars (Blackburne et al. 2006), different flux ratios in 
the lensed images were detected between the optical and 
X-rays. This dependency of the flux ratio on the energy 
has been interpreted as due to microlensing, thus justify- 
ing different amplifications as a function of the emitting 
region size (see also Jovanovic et al. 2008). On the basis 
of these evidences, we can consider that the lack of the 
gamma-ray echo may be due to different flux ratios be- 
tween the 2 components in radio and gamma-ray energy 
bands, i.e. it might be that one component (due to the 
flux anomaly) has negligible contribution in gamma-rays. 
Consequently intrinsic variation in combination with mi- 
crolensing may be also considered. 
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Fig. 1.— The 1-month E> 100 MeV light curve (2-day time-bin) of PKS 1830-211 (2010 October 8 - 2010 November 8). The two 
horizontal lines represent the 1-month and the 4-day 7-ray enhancements. On the right j/-axis we report the observed 7-ray SED. The 
downward triangles indicate the 2-cr 7-ray upper limits. The plot insert shows the soft X-ray light curve. 
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Fig. 2. — The observed SED of PKS 1830-211 for the 7-ray flare (red filled circles), for the 1-month enhancement (blue filled circles) and 
for the average state (black points). In detail, the black open circles are the non-simultaneous data from radio to hard X-rays (Pramesh 
Rao & Subrahmanyan 1998, Courbin et al. 2002, 7-year INTEGRAL survey, respectively). We report as open stars Planck data taken from 
the ERCSC of the first all-sky survey of Planck. The red downward arrows represent the near-IR SMARTS and the INTEGRAL/IBIS 
(20-40 keV) upper limits. Solid lines represent the SED models for the three 7-ray states (black for the average (EGRET), blue for the 
1-month, red for the flare). The models were magnified by a factor of ten for the lensing. 
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TABLE 1 

Model Parameters for 2010 7-RAY activity of 1830. 



population 


r 


B (Gauss) 


R(cm) 


K(cm- 3 ) 


7b 


7min 


ai 


a h 


8 


1st 


10 


0.7 


8 x 10 16 


100 


100 


35 


2.0 


2.6 


16 


2nd 


15 


0.4 


3 x 10 lb 


150 


500 


60 


2.0 


3.4 


20 



